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ABSTRACT 

We present a high-resolution millimetre interferometric image of the brightest SCUB A- 
selected galaxy from the Canada-UK deep SCUBA survey (CUDSS). We make a 
very clear detection at 1.3 mm, but fail to resolve any structure in the source. The 
interferometric position is within 1.5 arcsec of the SCUBA 850 /Ltm centroid, and also 
within 1.5 arcsec of a 44 /zJy radio source and a very faint, extremely red galaxy 
which we had previously identified as the submillimetre source. We also present new 
optical and infrared imaging, and infrared spectroscopy of this source. We model the 
overall spectral energy distribution and conclude that it lies within the redshift range 
2< z < 4.5. The submm/FIR luminosity of CUDSS14A is very weakly dependent on 
redshift within the constrained range, and is roughly 4xl0 12 L© (for Ho=75 and an 
assumed Arp220-like spectrum), which implies a star-formation rate ^1000 M yr . 
We derive an approximate gas mass of ~ 10 10 Mq which would imply the current star- 
forming activity cannot be sustained for longer than about 10 million years. With the 
present data however we are unable to rule out a significant AGN contribution to the 
total luminosity. 
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1 INTRODUCTION 

The commissioning of SCUBA on the JCMT has revolution- 
ized the field of observational cosmology and in particular 
the study of galaxy evolution, arguably making as big an 
impact as the release of the Hubble Deep Field did a few 
years earlier. 

Something of the order of 50 galaxies have now been 
discovered in the various deep SCUBA surveys (Smail et al 
1997, Hughes et al 1998 , Barger et al 1998, Eales et al 1999) 
and at least 30 percent of the far-infrared/submillimetre 
background discovered by COBE (Puget et al 1996, Fixsen 
et al 1998) has been resolved. Considerable effort is now 
being dedicated at ground-based telescopes in the radio, in- 
frared and optical, as well as Satellite observatories such 
as the Hubble Space Telescope,XMM-Newton and Chandra, 
to following up these sources, making identifications and ob- 
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taining redshifts, so that the history of star-formation in the 
Universe can be properly characterised free of the obscuring 
effect of dust. 

Accurate positions are of course essential for identifica- 
tions and follow-up observations, and the 15 arcsec SCUBA 
beam at 850 fim allows an uncertainty of several arcsec- 
onds even on a high signal-to-noise detection. Deep VLA 
imaging has been quite successful at providing identifica- 
tions and accurate positions of a number of SCUBA-selected 
sources (Smail et al 2000, Richards 1999). However, < 50% 
of SCUBA sources are detected with the VLA, even at ~ 
lO^iJy sensitivity levels, and in some cases even with a VLA 
detection there may be some ambiguity in the identifica- 
tion (Richards 1999, Downes et al 1999). The only way to 
truly confirm the identification is to make a higher resolution 
millimetre or submillimetre image using an interferometric 
array (Frayer et al 1998,1999,2000; Downes et al 1999). 

In this paper we present the results of an observa- 
tion with the IRAM interferometer on Plateau du Bure of 
CUDSS14A, the brightest source discovered in the CUDSS 
as reported by Eales et al (1999), along with new high res- 
olution J and K images obtained with the CFHT Adaptive 
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Table 1. Radio, Submm and ISO photometry 



Figure 1. The 1.3 mm map of CUDSS14A, contours are in steps 
of 0.5 mjy, the 1 sigma uncertainty, and extend from -2 crto 
+5cr. The cross marks the nominal SCUBA position on which the 
observation was centred. The effective beam shown is 2.7 X 2.0 
arcseconds with the major axis at a PA of 93 degrees. 



Optics system PUEO which we combine with archive HST 
images at 606 and 814 nm (roughly V and I), and also new 
infrared spectroscopy obtained on the United Kingdom In- 
frared Telescope (UKIRT). 



2 OBSERVATIONS 
2.1 IRAM Observations 

The IRAM data were taken on 25, 26 and 28 November 
1998. The array was in the 5D configuration, and the source 
was tracked for 8.4, 3.3 and 10.6 hours on each night re- 
spectively. Flux calibration was performed against CRL618, 
3C273 and MWC349 with 1637+574 and 1418+546 being 
used as phase calibrators. The observations were centred on 
the nominal SCUBA position of (J2000) 14 17 40.3 +52 29 
08 and observations were made at 1.3 mm (238.5 GHz) and 
2.8 mm (105 GHz) simultaneously. The data were correlated 
and reduced at the IRAM headquarters in Grenoble using 
the standard software. 

The 1.3 mm map is shown in Figure 1 (we note that 
Bertoldi et al 2000 have produced an essentially identical 
map from the same data). The object is very clearly detected 
at a level of 6 a, and the positional centroid is 14 17 40.37 
+52 29 06.8, or 0.6 arcsec in RA and 1.2 arcsec in declination 
away from the nominal SCUBA position. We estimate the 
total uncertainty of the IRAM position to be ± 0.3 arcsec 
(see Downes et al 1999 for a detailed discussion of the dif- 
ferent factors affecting the positional uncertainty of IRAM 
interferometer maps). There is no evidence for any exten- 
sion beyond the beam. We also examined the data for any 
evidence of spectral line emission but found none, however 
this does not place any meaningful limits on the redshift. 

There was no detection at 2.8 mm. The 1.3 mm flux 
and 2.8 mm upper limit are shown in Table 1, along with 
the SCUBA and ISO measurements. The overall flux dis- 
tribution, which is very similar to other SCUBA-selected 
galaxies (e.g. Ivison et al 1999), is plotted in Figure 2. 



Wavelength 


Flux (mjy) 


reference 


6 cm 


0.044 ±0.004 


Fomalont et al.1991 
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this paper 
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2.94 ±0.49 


this paper 
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Eales et al 1999 
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Figure 2. The overall radio through optical observed flux distri- 
bution of CUDSS14A. 



2.2 New CFHT Observations 

Observations were made with the PUEO Adative Optics sys- 
tem on CFHT during the nights of April 27-28 1999, using 
an R=15.3 guide star 26 arcsec WNW of the target. The 
detector was a 1024 2 HgCdTe array sampling the image at 
0.035 arcsec/pixel. Secure detections were obtained in a to- 
tal of 150 minutes integration in J and 100 in K'. 

The new infrared images have been co-registered with 
archival HST images in F606W and F814W (the source is in 
the so-called 'Westphal' survey field), resampling to a pixel 
scale of 0.1 arcsec/pixel. A montage of the F606W, F814W, 
J and K' images is shown in Figure 3. The morphology of the 
source appears compact and indistinct. The FWHM of the 
images are about 0.4 arcsec in all of the bands, substantially 
larger than the PSF in the HST data. The PSF is poorly 
constrained in the Adaptive Optics images but is most likely 
substantially smaller than 0.4 arcsec. 



Table 2. Multi-aperture optical/infrared photometry (AB mags) 



Filter 


1 arcsec 


2 arcsec 


4 arcsec 


F606W(HST) 


27.132±0.314 


27.011±0.351 


26.063±0.263 


F814W(HST) 


25.228±0.127 


25.023±0.128 


24.574±0.128 


J(CFHT) 


23.230±0.094 


22.866±0.124 


22.686±0.215 


K'(CFHT) 


21.423±0.047 


21.195±0.072 


21.495±0.197 
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Figure 3. Images at F606W, F814W, J and K' clockwise from top 
left respectively of CUDSS14A. The images are all co-centred on 
the extremely red optical source, and are all 10 X 10 arcseconds. 

The brightness of CUDSS-14A in F606W, F814W J and 
K' in apertures of 1, 2 and 4 arcsec diameter are listed in 
Table 2. The colours are broadly consistent in the different 
apertures but with a trend to be bluer in larger apertures. 
This new photometry agrees well with earlier measurements 
in larger apertures by Lilly et al (1999). Optical spectroscopy 
was also attempted but no detection was made. 

2.3 UKIRT Observations 

CUDSS14A was observed with the CGS4 spectrometer on 
the nights of April 13 and 16 1999 using the 300 mm focal 
length camera and 40 lines/mm grating. The 2.4 arcsec slit 
was centred on the radio position, and the observing and 
data reduction procedures described by Eales and Rawlings 
(1993) were used. In the K and H bands we observed in first 
order (X/SX = 400) for 160 minutes, and in J we observed 
in second order (X/SX — 800) for 80 minutes. 

Continuum was detected at K, but not at J or H. No 
spectral lines were seen in any of the bands and we esti- 
mate the upper limits on the lines as 1.1 x 10~ 19 W m -2 
in the region 2.02-2.3 /im, 1.3 x 10~ 19 W m~ 2 in the region 
1.5-1.8 /im, and 2.8 x 10 -19 W m -2 in the region 1.1-1.3 
/jm. To derive these limits we have assumed that a spectral 
line would have a spectral width equivalent to a velocity 
dispersion of 200 km s . These limits are slightly more 
conservative than 3<r limits, since we have enough spectral 
range that we might expect a 3u fluctuation by chance. The 
limits have been chosen to be slightly larger than the largest 
feature in each spectral range. 



3 DISCUSSION 

3.1 The identification 

The IRAM position is 1.4 arcsec from the nominal SCUBA 
position and is also only 1.2 arcsec from the 5 GHz position 



of Fomalont et al (1991) and 1.5 arcsec from the optical ID 
(Lilly et al 1999). Since there is no other candidate in each 
band than the one we show and all the positions lie easily 
within their respective 2a astrometric uncertainties of each 
other (in fact the IRAM, radio and optical positions are only 
just outside their respective 1-a uncertainties) we believe the 
identification to be extremely robust. 



3.2 The Redshift 

In Lilly et al (1999) we suggested a redshift in the range 
z=2-3, based on the rather uncertain optical and infrared 
colours available at that time. We can now use our new 
data to provide further photometric constraints. 

3.2.1 Radio and Submillimetre 

Carilli and Yun (1999; hereafter CY99) suggested that the 
strong radio-FIR correlation could be used to provide con- 
straints on the redshifts of SCUBA survey sources. The orig- 
inal CY99 relationship is based on a model, but in a follow- 
up paper (Carilli and Yun 2000, hereafter CY00) they used 
observational data from 17 galaxies observed in the submil- 
limetre by Lisenfield, Isaak and Hills (1999) to produce an 
empirical relationship, which tends to result in somewhat 
lower redshifts than their original one. Both these results 
have now been superseded however by Dunne et al (2000) 
who have made a SCUBA survey of 104 nearby galaxies. 
Using Dunne et al's empirical relationship between 850 fim 
and 1.4 GHz flux density, including the intrinsic scatter to 
derive an uncertainty, the measured 1.3 mm to 850 /im spec- 
tral index of 3.2 and the 1.4 to 5 GHz radio spectral index 
of -0.2±0.3, we obtain z~2.2±0.5 for CUDSS14A. For com- 
pleteness we note that CY00 would result in z~2.9±0.7 and 
CY99 in z~5.0±1.0. The rough consistency of the CY00 and 
Dunne et al results and the inconsistency of the CY99 pre- 
diction demonstrates both the uncertainty in this technique 
(see also Blain 1999) but also the manner in which it is 
rapidly improving. 

We can also use the fact that we have upper limits at 450 

and 2.8 mm to provide further constraints. To do this 
we take two extreme template spectra, namely the nearby 
starburst galaxy M82, which has a dust temperature of 48K 
and is optically thin throughout the FIR and submm region 
(Hughes et al 1994), and Arp220, which has a fitted dust 
temperature of 62K and is optically thick shortward of 200 
fiva (see e.g. Downes, Solomon & Radford 1993). 

As demonstrated in Figure 4, CUDSS14A cannot be at 
low redshift, for either template, otherwise we would have 
detected it at 450 /im. Conversely it is unlikely to be at 
very high redshift, for either template, otherwise we would 
have detected it at 2.8 mm. In fact the detections and upper 
limits constrain the redshift to the range 2< z < 5. 

It would be possible to allow a lower redshift if the emit- 
ting dust were much colder. In order to be consistent with 
the 1.3 mm detection and 450 fim upper limit however, for 
z<l the temperature would have to be < 20K. Although 
this would reduce the source luminosity somewhat (depend- 
ing on the unobserved FIR spectral shape), it also increases 
the required gas and dust mass to > 1O 11 M0, which is rather 
excessive. Interestingly, the 450/1300 fim flux constraint 
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Figure 4. Comparison of the mm/submm spectrum of CUDSS14A with that of (a) M82 and (b) Arp220 at z=2.0 (dots), 2.5 (solid) and 
5 (dashed) with arbitrary shifts in absolute flux level. Note that the M82 spectrum is not consistent with the data for z=2.0 but is just 
consistent for z=5, whereas Arp220 is just consistent at z=2.0 but is not consistent for z=5. Thus for these templates we can constrain 
the redshift to the range 2.0 to 5.0. 



for CUDSS14A is identical to that found for HDF850.1 by 
Downcs ct al (1999) and they reach the same conclusion, 
that the source is unlikely to be at such low redshift. We 
also note that no sources as cold as 20K are known at low 
redshift and furthermore, that higher luminosity sources also 
tend to be warmer. Even at the lower end of any possible 
luminosity range it seems very unlikely to us that such a 
luminous object would be very cold. 

3.2.2 Optical and Infrared 

CUDSS14A is an extremely red object (ERO) with (V - 
K)ab ~ 6 corresponding to a Vega-normalized colour of 
(V-K) ~ 8. Other EROs have been identified with SCUBA 
sources (Smail et al 1999), and Barger, Cowie and Richards 
(2000) have shown that many very red radio source identi- 
fications are detectable with SCUBA. An interesting conse- 
quence of the very red colours of a source like CUDSS14A 
is that such galaxies would rapidly become extremely faint 
if placed at higher redshifts (Dey et al 1999). This is illus- 
trated in Figure 5 where we show the expected V, I, J and 
K magnitudes as a function of redshift assuming that the 
real source is (a) at z—2 and (b) at z=4. In the former case, 
the source would become very hard to detect at wavelengths 
below 1/im at redshifts as low as 3. The implications of this 
are that SCUBA sources for which no optical identification 
can be found even to extremely faint levels such as reached 
by the Hubble Deep Field are not neccessarily at very high 
redshifts, z^>3. 

The UK colours of CUDSS14A are consistent with an 
Sb galaxy at z=2.5 (as suggested in paper II), however the 
VIJ colours are also consistent with an Sa galaxy at z=1.5, 
and the overall spectrum is so extremely red it is not very 
different from a power law. Therefore we cannot make any 
firm conclusions based on the colours, but most probably z 
is less than 4.5 based on the V detection. 




2 4 S 8 10 13 
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Figure 5. The expected variation of the V, I, J and K magnitudes 
(in AB) with redshift of a source identical to CUDSS14A derived 
assuming that CUDSS14A as observed lies at z = 2 (solid lines 
starting at the solid symbols at z = 2) and at z = 4 (dashed lines 
starting at the open symbols at z = 4, no V curve is plotted in 
this case since the Lyman break should suppress V very rapidly 
above z ~4). The horizontal dashed area in grey represents the 
approximate magnitude limits of the HDF. The diagonal grey 
lines represent more typical limits to ground-based imaging and 
are shown for the two sets of magnitudc-rcdshift predictions. An 
f! = 1 Universe has been assumed. 

If CUDSS14A were at z ~2.2-2.5 then Ha would appear 
in the range 2.1 to 2.3 /im, and the [OIII] 5007 line would 
appear at 1.60 to 1.75 fim. No lines are observed. However 
we note that spectroscopic observations of nearby ULIRGs 
by Kim, Veilleux and Sanders (1998) imply the lines are 
suppressed compared to the level expected from the FIR- 
derived star-formation rates in these objects, and in fact if 
the ratio of optical/infrared line emission to FIR luminos- 
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ity in SCUBA-selected galaxies is similar to that found in 
nearby ULIRGs then 8-10 m telescopes will be required to 
detect them. Therefore the non-detection of any lines cannot 
be used to rule out a redshift of 2-3. 

The optical/infrared, radio and mm/submm data there- 
fore all argue strongly for a redshift 2.0 < z < 4.5 with some 
preference for a redshift around 2.2 to 2.5. 



~10 Mq. The uncertainty in this estimate is probably at 
least a factor 5 in each direction however, given the un- 
certainty in dust temperature, optical depth, grain proper- 
ties, redshift and cosmology. A galaxy of this mass could 
only sustain the derived star formation rate for a period of 
around 10 million years. So these results are all consistent 
with CUDSS14A being a gas and dust-rich galaxy undergo- 
ing a short-lived period of extensive star-formation. 



3.3 AGN Heating ? 



A source at redshift z>2 with a dust temperature around 
50K has a minimum possible source size of 0.1 arcsec to 
produce an observed flux of 2.9 mjy at 240 GHz, which 
at z~2-5 corresponds to roughly 1—2 kiloparsec, which is 
the typical size for circumnuclear starburst regions in lo- 
cal ULIRGs. Another way of looking at this is to say that 
if the emitting region in CUDSS14A is a typical circumnu- 
clear star-forming region and it is at z~2-5, then the submm 
emission must be quite close to being optically thick. 

Applying the same argument to the question of whether 
the FIR luminosity could be produced by AGN heating, we 
can say immediately that the submm-emitting dust cannot 
predominantly be heated directly by a central AGN, since in 
order to be at a temperature of 50K or below, for standard 
dust properties (e.g. Hildebrand 1983) the grains would have 
to be ~ 5 kpc from a central AGN of luminosity 4xl0 12 Lq 
and would therefore have been resolvable. Conversely if opti- 
cally thin dust were only 1 kpc from such a central source its 
temperature would be ~ 200K which is not consistent with 
the ISO and 450 /im non-detections. However we cannot rule 
out heating via radiative transfer through a large column of 
dust, with the observed mm/submm emission arising only 
close to the outer surface. 

The launch of the Chandra and XMM X-ray satellites 
offer the opportunity to observe AGN emission directly in 
these sources if it is present and should help to solve this 
question definitively. Initial Chandra observations presented 
by Mushotsky et al (2000) suggest that the hard x-ray back- 
ground has been resolved into heavily obscured active galax- 
ies, for which the SCUBA-selected objects are a plausible 
counterpart. However the limited data so far available from 
Chandra on SCUBA sources (Fabian et al 2000) suggest 
that these are more like starburst galaxies, unless they are 
Compton-thick and the scattered X-rays are weak or also 
absorbed. 



3.4 Source properties 

We can derive a total luminosity for CUDSS14A of ~3- 
6xl0 12 Lq for 2<z<4.5. This luminosity estimate is for an 
Arp220-hke spectrum where the total luminosity is com- 
pletely dominated by the Submm/FIR. For an M82-like 
spectrum the total luminosity could in fact be considerably 
higher. For a standard initial mass function, this luminosity 
corresponds to an ongoing star-formation rate of ~1000 M© 
yr -1 , which is hard to sustain for a cosmologically signif- 
icant time, unless the IMF is very heavily biased towards 
massive stars. 

We can also estimate a total gas mass using the tech- 
niques laid out in Hildebrand (1983) and Gear (1989) (see 
also Hughe, Dunlop and Rawlings 1997) to give roughly 
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